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Abstract 

We present calculations of rotational absorption spectra of the molecules HCN and DCN in 
superfluid helium-4, using a combination of the Diffusion Monte Carlo method for ground state 
properties and an analytic many-body method (Correlated Basis Function theory) for the excited 
states. Our results agree with the experimentally determined effective moment of inertia which 
has been obtained from the J = — > 1 spectral transition. The correlated basis function analysis 
shows that, unlike heavy rotors such as OCS, the J = 2 and higher rotational excitations of HCN 
and DCN have high enough energy to strongly couple to rotons, leading to large shifts of the 
lines and accordingly to anomalous large spectroscopic distortion constants, to the emergence of 
roton-maxon bands, and to secondary peaks in the absorption spectra for J = 2 and J = 3. 

PACS numbers: 05.10.Ln, 05.30.Jp, 33.20.Bx, 33.20.Ea, 34.30. +h, 36.40.-c, 67.40.Yv 
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I. INTRODUCTION 



In microwave helium nanodroplet isolation spectroscopy experiments, Conjusteau 
et. al. Q have measured the rotational excitation energy J = ^ 1 of HCN and DCN 
embedded in ^He clusters. Their results show a reduction of this excitation energy by fac- 
tors of 0.815 and 0.827 with respect to gas phase HCN and DCN, respectively. Infrared 
spectroscopy experiments of HCN by Nauta et. al. 2] yield similar results from analysis of 
the ro- vibrational excitation of the C-H stretching mode, namely a reduction of 0.795 in the 
J = — ^> 1 energy. These fractional reductions are considerably smaller than those observed 
for heavier molecules such as SFg and OCS, where reductions by factors of ~ 0.3 are seen 
The gas phase rotational constants, B = 1.478222 cm~^ for HCN and B = 1.207780 cm~^ 
for DCN, are also much larger than the corresponding values for the heavier molecules (e.g. 
B = 0.2029 cm~^ and 0.0911 cm~^ for OCS and SFg respectively). The widely observed re- 
duction in B is understood to be due to the interaction of the molecule with the surrounding 
^He awn. Q. Fo. the heav.. .„,ecu,e= . ha= been found that calcnlat.on. based on the 
microscopic 2-fluid theory can reproduce the effective rotational constant B^s 0,0]. For 
some heavy linear rotors, a semiclassical hydrodynamical analysis that combines a classical 
treatment of the molecular rotation with a quantum calculation of helium solvation density 
approximately reproduced the moment of inertia increase measured in experiments (see ta- 
ble I in Ref. although no agreement is found for the octahedral SFg molecule 0,0,0]. The 
hydrodynamic contribution to the effective moment of inertia is found to be considerably 
decreased when the molecular rotation is treated quantum mechanically |lo| . 

These models for heavier molecules are based on analysis of partial or complete adiabatic 
following of the molecular rotational motion by helium and cannot describe the dynamics of 
light rotors like HCN and DCN in helium for which adiabatic following does not hold jlll ]. 
Furthermore, infrared spectra of HCN [2] and acetylene, C2H2 jl2|], and other light molecules 
show a small splitting of the ro- vibrational R(0)-line which cannot be accounted for by these 
theoretical approaches that focus on the calculation of the helium-induced increase of the 
moment of inertia. In Ref. |l3|, a detailed investigation of the effects of the finite ^He envi- 
ronment on rotational excitations showed that both hydrodynamic coupling of translational 
and rotational motion and the anisotropy of the effective interaction between the molecule 
and the finite ^He cluster, can result in splitting of the -R(O) spectral line (corresponding to 
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the J = — s> 1 transition) into M = and M = ±1 contributions. However, the observed 
hne shapes could not be explained in the case of HCN, although good agreement was found 
for the lineshape of the R(0) transition of the heavier rotor OCS. 

The molecules HCN and DCN are light rotors, possessing large zero point motion. There- 
fore, calculation of the ground state already requires a full quantum mechanical treatment 
of the molecular rotations 14]. Furthermore, for the rotational excitations, the spacing 
between the rotational energy levels is large, of similar magnitude as the roton energy of 
bulk helium (the roton gap is 8.7K This introduces the possibility of direct coupling 

between the roton states and the molecular rotational levels of light molecules. The cou- 
pling between phonons in ^He and molecular rotational levels was analyzed perturbatively 



in Ref. |l6|, where it was shown that the lower density of phonon states in ^He relative to that 
of particle-hole states in ^He leads to a much lower coupling of molecular rotational transi- 
tions to excitations of the quantum liquid for the Bose system. This provided a rationale 
for the observation of sharp rotational lines in infrared molecular spectra in the bosonic ^He 
environment, but not in the fermionic ^He environment. The specifics of the dispersion rela- 
tion in ^He were not incorporated in this perturbative analysis. In particular, only a linear 
phonon spectrum was employed, and the maxon and roton excitations were not taken into 
account. To allow for the possibility of coupling to maxons and rotons when calculating the 
response to the motion of these light molecules, it is evident that helium cannot simply be 
treated as a classical frictionless fluid possessing long wavelength hydrodynamic modes, nor 
by a quantum fluid possessing quantized phonons with linear dispersion. We must therefore 
describe the coupled dynamics of the molecule and the strongly correlated helium quantum 
fluid with true quantum many-body theory, i.e., in principle we must solve the + 1 body 
Schrodinger equation. 

Quantum Monte Carlo is one such fully quantum approach. Zero temperature Quantum 
Monte Carlo calculations for the rotational motion of a molecule doped into a finite cluster of 
Keliu. have been ean^ed out successfully fo. a va*.y of s,nall molecules flU QQ. In 

which excitation energies have been obtained with the POITSE (Projection operator imagi- 
nary time spectral evolution) method [19] or similar spectral evolution approaches For 
HCN in helium clusters, rotational constants -Bgfr have been obtained for clusters consisting 
of up to = 25 ^He atoms Q]. However, in marked contrast to the heavier OCS and SFg 
molecules for which the large droplet value is arrived at well before the first solvation shell 
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(A^ ~ 20) is complete |12I,I llSl |20|, convergence to the experimental value of B for HCN in 



large ^He clusters of several thousands of ^He atoms was not yet found at = 25. This very 
different behavior accentuates the distinction between a light and a heavy rotor, and suggests 
that different physics may underlie the reduction in B for a light molecule. This, combined 
with the technical difficulties of making POITSE calculations for large clusters and hence to 
follow the small cluster B value to convergence with increasing size, motivates development 
of a more analytic approach that is suitable for implementation in the bulk limit. With 
an analytical approach, implementation is generally simpler in the bulk limit N ^ oo than 
in a large finite system because of the higher symmetry. However, it is then necessary to 
recognize that derivation of an analytical model from the + 1-body Schrodinger equation 
is necessarily approximative and this may affect the quantitative accuracy of our results. 

The method we apply here to the rotational dynamics of HCN in helium (1 linear molecule 
+ spherical atoms) is a combination of the Correlated Basis Function (CBF) theory (also 
called (time-dependent) hyper-netted chain / Euler-Lagrange (HNC/EL) method) with dif- 
fusion Monte Carlo (DMC) ground state calculations. The CBF method can be formulated 
as an energy functional approach to solving the many-body Schrodinger equation. In con- 
trast to the formulation of Density Functional Theory (DFT) that is generally applied to 
helium systems 0, CBF theory eliminates the need for a semi-empirical correlation 
energy functional by expressing the energy functional not only in terms of the one-body 
density, but also in terms of pair-densities and, if necessary for quantitative agreement, also 
of triplet-densities. Similarly to DFT, the stationary version of CBF yields the ground state 
energy and structure, while the time-dependent extension of CBF yields excited states. 

As an analytic approach to the many-body problem, CBF theory requires relatively little 
computational effort to solve the equations of motion, once these have been derived. The 
CBF method yields ground state quantities such as the ground state energy, the chemical 
potentials, and the pair distribution functions. Calculation of excitations in CBF yields 
not only excitation energies, but also the density-density response function, and from that 
the dynamic structure function for pure '^He and the absorption spectrum of the dopant 
molecule, as will be shown explicitly below. Although CBF theory is not an exact method, 

evant to ^He systems. 



quantitative agreement has been found for a variety of quantities re 
These include ground state and collective excitations in bulk ^He jsS, 

films 2^ 21, 2^ and clusters [3, ^ImI, and translational motion of atomic impurities 



2J], excitations in "^He 
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in ^He j32|. This generates confidence that CBF theory may allow quantitative calculations 
of the rotational dynamics of dopant molecules for rotational energies in the range of the 
phonon-maxon-roton regime. In the case of HCN/DCN, this means for quantum numbers 
up to J = 3. Combining CBF theory for excitations with exact ground state quantities 
calculated by DMC may also be expected to improve the accuracy of CBF excitations. 

In this paper, we restrict ourselves to the simpler case of HCN/DCN in bulk "^He, where 
translational symmetry is preserved and the analytic CBF calculations become correspond- 
ingly simpler, as opposed to HCN/DCN in ^He clusters. This therefore precludes the calcu- 
lation of inhomogeneous line-broadening and possible line splitting caused by the inhomo- 
geneous environment of the cluster The structure of the paper is as follows. 

The derivation of the rotation excitation spectrum of a single linear molecule in bulk 
helium is presented in section [H] This analysis is related to the derivation of the translational 
excitation spectrum of an atom coupling to the phonon-roton excitations in bulk helium, 
which has been discussed in detail in Ref. 0- In sections III Al and IIIB[ we derive the 
CBF equations for the excitation spectrum of a linear rotor solvated in bulk ^He and for 
its absorption spectrum, respectively. In section III C\ we describe how the CBF theory is 
combined with Diffusion Monte Carlo (DMC) calculations of ground state input quantities 
in order to obtain excitations and absorption spectra. 

In section IIIII we report the DMC results for the ground state quantities for HCN and 
DCN. Section IIVI describes the results of the "marriage" of CBF and DMC for rotational 
excitations. We present and discuss the excitation energies in section HV Al and the absorp- 
tion spectra in section IIVBI In a direct analogy with the definition of an effective mass of 
an effective free particle from the momentum dependence of the excitation energy E{q), it 
is common to obtain an effective rotational constant i?efr of an effective linear rotor from 
analysis of one or more spectral transition energies. This was done in Ref. Q using the ex- 
perimentally observed R(0) (J = ^ 1) spectral line. We present our corresponding results 
for Bcs obtained from the calculated J = —>■ 1 transition energies in section II V (^1 Our 
calculated values of Bes are in good agreement with the experimental values, indicating that 
the reduction in B^s relative to the gas phase B for these light molecules derives primarily 
from coupling to the collective modes of ^He. This is a very different situation from that 
for heavier molecules, where the reduction in B^s derives from coupling to some local he- 
U.n density .ha* ad,aU«ca..y follow. t.e .o.ecu.a„ota«o„Q , a phenomenon that may 
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be formally regarded as coupling to ^He modes which are localized around the molecule. 
The present analysis thus indicates that there is indeed a different physics responsible for 
the reduction in rotational constants for light molecules than for heavy molecules in ^He. 
A second significant feature of the CBF results is that, although we find that the energy 
spectrum still has the same symmetry as a linear rotor, i.e., there is no splitting of the M 
states within a given level in bulk ^He, it is evident that nevertheless the J dependence of 
the rotational energy E{J) deviates considerably from that of an effective rigid linear rotor 
spectrum, B^ffJ^J + l), when a fit of i?efr to more than one J level is made. In particular, we 
find that the coupling to the roton and maxon collective excitations for higher J levels gives 
rise to extremely large effective "centrifugal distortion" terms that modify this rigid rotor 
form. The analysis of this deviation from the rigid rotor spectrum is discussed in detail in 
section IIVDI 

Lastly, in section IIV Ft we introduce a pseudo-hydrodynamic model that includes only 
phonon modes of ^He but no maxon/roton modes in the CBF calculation. This provides a 
reference point that allows us to independently assess the effect of the maxon/roton exci- 
tations on molecule rotations. The changes in the effective rotational constants B and D 
relative to the gas phase values deriving from this pseudo-hydrodynamic dispersion model 
are much reduced relative to the corresponding changes found with the true disperson curve 
for ^He, and the value of B^s is no longer in such good agreement with the experimen- 
tally measured value {Dcs was not experimentally accessible for HCN in the experiments 
to date 

HQ). 

This provides additional evidence for the critical role of the maxon/roton 
excitations in the reduction of B^s for HCN. We summarize and provide conclusions in 
section M 



II. THEORY 



The CBF method is a microscopic quantum theory for the ground state and excitations 
of a many-body system. By "microscopic" we mean that there is no input other than the 
Hamiltonian, and the output quantities are expectation values with respect to the ground 
state or an excited state, such as energy, density etc. In practice, approximations are neces- 
sary in order to render the CBF equations soluble. We will point out these approximations 
as we introduce them. 
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In our case the Hamiltonian for '^He atoms with coordinates rj, i = 1, . . . ,N and a 
hnear molecule at position fq and orientation Q = {6, 0) in the laboratory frame takes the 
form 

H = BL'- —Vl + Y: Vxiro - r„ Q) + , (2.1) 

i=l 

where B is the rotational constant of the free linear rotor, L is the angular momentum 
operator, M is the mass of the rotor, and Vx is the molecule-^He interaction potential. For 
HCN-He, we use the 1E8 potential of Atkins and Hutson js^ obtained from fitting to ab initio 



calculations of Drucker et al. 



For DCN-He, we use the same potential (same equilibrium 

I) and merely transform the Jacobi 



nuclear positions, rcH = 1.064A and rcN = 1.156A 
coordinates (r, a) to take into account the change of the center of mass, r = |ro — rj| is the 
helium distance from the molecule center of mass, a is the angle between the vector r and 
the molecular axis, measured from the hydrogen end of the molecule. 
The operator Hb is the pure helium Hamiltonian 

= -— E V? + E ^He(|r. - r, I) , (2.2) 

i=i i<j 

where m is the mass of an ^He atom and l^He is the ^He-^He interaction, for which we use 
the potential of Ref. iSf 



The CBF method has been explained in detail in a number of papers 



H H Q, 



4ll |. therefore we limit ourselves to giving only a very brief overview here. The starting point 



is to obtain the ground state wave function of the A^ + 1-body system, here 

^0 = ^o(ro,ri,...,r7v,l^). (2.3) 

In the framework of CBF theory, \l/o is expressed in a Jastrow-Feenberg form, i. e. expressed 
in terms of correlations: 



^0 = exp^ 



J2u2iri,rj)+ E M3(r„rj,rfc) + . . . 

i<j i<j<k 
N -| 

+ E«f(ro,r„fi) +E«f(ro,ri,rj,(]) + ... . (2.4) 

i=l i<j 

Here the molecule is referred to as X, with center of mass translation coordinate Tq and 
angular orientation Q defined above. The definition of the ra-particle correlations u„ and 
is made unique by requiring that Un vanishes if one of its coordinates is separated from 



the rest. Furthermore, it provides an exact representation of the ground state when all 
correlations up to = are summed, i.e., up to un- However, even for the strongly 
correlated ^He ground state, correlations between up to just three particles are sufficient to 
obtain quantitative agreement of the energy and the pair distribution function 38] of bulk 
■^He with experiments and Monte Carlo simulations. The correlations m„ are obtained by 
solving the Euler-Lagrange equations, which can be written formally as 

5E , , 

r = (2.5) 

bUrSji, . . . ,r„,) 

where E is the expectation value of the Hamiltonian, (\l'o|i/|\l'o), and < 3. The resulting 
Euler-Lagrange equations (j2.5j) are coupled non-linear integro-differential equations and can 
be solved iteratively. Derivation of a formulation of equations ()2.5|) that is appropriate for 
numerical solution can be found in Ref. |40|. However, in the present "marriage" of CBF and 
DMC, solution of eqns. ()2.5p is not necessary since the ground state properties are calculated 
by DMC. 

A. Excited states 

The primary aim of this article is to employ CBF theory in the search for excitations of 
the molecule-helium system. The excitations can be obtained by generalizing the ground 
state form, Eq. ()2.4|1 . to time-dependent correlations, i.e. ui = M2(ri, ^2] t), etc. By allowing 
a time-dependent external perturbation potential V^'^^^\rQ,^l;t) to act on the molecule, we 
can then use linear response energies involving motions of 

the molecule. Linear response relies on the knowledge of the ground state, which we assume 
to have calculated according to the above recipe or by other means (see section III C|) and 
which is then weakly perturbed. The perturbed wave function can therefore be written as 

g5C/(i)/2^^ 

where the excitation operator 6U (t) is given by 

N 

5U{t) = 6ui{ro, n-t) + J2 ^Mt^o, fl; t) . (2.7) 

1=1 

Note that we have dropped the X superscript from the two-particle molecule-helium cor- 
relation U2. We will continue to do this from here on, using the presence of the molecular 
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coordinates Tq and Q in Sui^tq, Q; t) and 5u2{vq, Vi, Vt; t) to distinguish helium-molecule from 
helium- helium correlation terms. Note also that unlike the ground state wave function 
the excited state '^{t) does not possess the translational and rotational symmetry of the full 
Hamiltonian H . Therefore it is convenient to separate 6U{t) into a one-body term 6ui and 
2-body correlations 6u2- Time-dependent correlations between more than 2 particles can 
also be formally written down and added to Eq. ()2.7p . but would give rise to numerically 
intractable equations. Consequently, we restrict ourselves to two-body correlations in 6U{t) 
here. 

The two terms in Eq. ()2.7j) give rise to two Euler-Lagrange equations that are obtained 
by functional minimization of the action integral 

r*2 , d 



6C = 6 [' dt + -i^— =0 



(2.8) 



with respect to 5ui and 5u2. The action integral of an __s5herical impurity (^He and atomic 
hydrogen) in bulk ^He can be found in Refs. 



dbd and |42|. We shall refer to this reference 
integral for a spherical impurity as Cq. For a linear molecule in helium, the situation is 
complicated (i) by the additional rotational kinetic energy term in the Hamiltonian, 

1 d , . B ^ 1 



rrrot 
-^0 



BL' = -B I ^A(sin^^) + 



(2.9) 



and (m) by the breaking of the rotational symmetry of the ground state distribution of ^He 

1 \ 

atoms around the molecule. Similarly to the derivation of Cq [323, ^'^'^ expansion 



of C to second order in 5U: 

rt2 



c = c 







-B I dt< 

4 iti 



X 



|ae(5ni(0,fi)|2 + 



2n 



sin'^ 



dMQdl p2(0, 1,1^) 



(ae5ut(0,fi))(9e5u2(0,l,r])) + c.c 



+ 



+ 



sin0 j \ 
d0dndld2 p3(0, 1,2,1]) 



sin 6 



c.c+\de5u2i0,l,^)\' 



\d^Su2{0,l,n)\ 



2-1 



sin' 9 



{de6u;{0, 1, n)) {deSu2{0, 2, Q)) + c.c 



+ 



sin 6 



sin 6 



t2 



dt / dOdQ pi[6U]{0,Q)V^'''''\0,Q) 



where for simplicity we abbreviated the functional arguments by i, and have omitted the 
time argument. P2(0, 1, ^) and P3(0, 1, 2, fl) are the ground state probability densities of one 



and two ^He atoms around the molecule, respectively, defined as 

P2(0,l,fi) = -^y c/2...c?iV 1^0(0,1,..., iV,f])|2 (2.10) 

P2(0,l,2,fi) = ^^^i^^^ 1^3... rfiV 1^0(0, l,...,iV,fi)|2, (2.11) 

where M is the normalization integral of \E'o. 

pi [5U] (0, r]) = + 3fJe5pi (0, n) (2.12) 

is the time- dependent probability density of the molecule expanded to first order in 5U, 
where we have defined the complex density fluctuation 

(5pi(o, n) = p^Suiio, n) + J di p2(o, 1, n)5u2{o, 1, n) , (2.13) 

and p^ = 1/V and p = N/V are the constant ground state densities of the molecule and of 
the ^He atoms, respectively, in the normalization volume V. 

The Euler-Lagrange equations 6C, i.e. the 1-body and 2-body equations 

0, . ./^ , =0. (2.14) 



describe the time-dependent response of the system to the perturbation \/(^^*). The linear 
response is obtained by linearizing the equations in terms of the corresponding correlation 
fluctuations 6ui and 6u2. In the following, we bring these equations into a form where 
the time-dependent density fluctuation 5pi(0,f2) is a linear functional of V^°^'^\ From that, 
excitations are derived by setting = 0. 

We first eliminate 6ui{0,Q) in favor of the (complex) one-body density fluctuation 
eq. (j2.13|) . Then the linearized 1-body equation of motion can be written 

it']- ^) + ^) + "^ofiO, n) - 2\/(-*)(0, fi) = , (2.15) 

sm6' \d^J 

where weak V^'^^^^ is any perturbation acting only on the molecular degrees of freedom and 



j is the translational current fluctuation. This is deflned in Ref. |32 and need not concern us 
for rotational excitations, as we will see explicitly below. In analogy to j"^, we have deflned 
the "rotational" current fluctuation 

r(0,fi) = (""'f^'] 5pi(0,fi) - fd2 6u2{0,2,n) hf^'] P2(0,2,fi) . (2.16) 

10 



Spi{0,Q) is the time-dependent density fluctuation deflned in eq. ()2.13|) . while, due to the 
hnearization, p2{0,2,Q) is the ground state pair density. The second and third terms of 
eq. (j2.15j) stem from the variation 6Co/6ul{0,Q). The density fluctuation 6pi couples via j*" 
and j"^ to the 2-body equation in eq. (j2.14|) . The 2-body equation is more lengthy: 

= ^9,sin^p,(0,l,l])9,^^l%^ + -^9,p,(0,l,(])a,^M^ (2.17) 
sm sm 9 p^ 

+ -^de sin^p2(0, 1, n)de5u2{0, 1, n) + -^-rS^PslO, 1, n)d^5u2{0, 1, n) 
sm 6 sm 6 

+ fd2 —-desmep3{0,l,2,Q)de6u2{0,l,Q)+ f d2 —^d^p3{0,l,2,n)d^6u2{0,l,n) 
J smd J sm (J 

- [d2 ^-^desmep2{0,l,^)deP2{0,2,n)Su2iO,l,n) 
J p^ smt7 

- fd2 ^-L^d^p2{o, 1, n)d^p2{o, 2, n)5u2{o, i, n) 

J p^ sm 6 

r /'trans 

+ ^M^2(0,l,fi) + Vi-J2(0,l,fi) - ° - 2p2(0,l,fi)\/(^^*)(0,l]). 

du2{0, 

SC^q'^^^ /6u2{0,l,Q) represents the terms related to the translational degrees of freedom of 



the molecule, the derivation of which again can be found in Ref. |3j. The (complex) 2-body 



density fluctuation 5p2(0, 1, Q) can be expressed as a functional of 5pi(0, Q) and 5^2(0, 

Sp2{0, 1, Q) = p2(0, 1, Q)^-M^ + p2(0, 1, Q)6u2i0, 1, ft) (2.18) 

P 

+ j d2 (p3(0, 1, 2, fi) - -^p2(0, 1, fi)p2(0, 2, fi))(5M2(0, 2, fi) , 

which follows from the deflnition ()2.10j) . The ^He-current fluctuation J2 induced by the 
rotating molecule is deflned as 

J2(0, 1, 9) = -^P2(0, 1, n)Vi5n2(0, 1, fi) . (2.19) 
2m 

Unfortunately, solution of the coupled set of equations ()2.15p and ()2.17p is not feasible 
without further approximations to the 2-body equation ()2.17|1 . In a rather drastic approxi- 
mation step, we therefore expand the pair distribution 

S(0.1.n) = ^^^%i^ (2.20) 

p p 

about unity. This is commonly referred to as the "uniform limit" approximation IstJ. It 
has the advantage of leading to a particularly simple expression for the excitation energies 
in terms of a self energy correction, and has been used in many CBF calculations of excited 
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states of ^He and impurities in ^He. Therefore, this uniform hmit approximation is our 
first candidate for simphfying eq. ()2.17p . We discuss the extent of the vahdity of this 
approximation in section ITTTl where we present our DMC result for g. 

When apphed to the equations of motions (j2.15|) and (j2.17|) . the uniform hmit approxi- 
mation amounts to replacing the pair distribution function (7(0, 1, Vt) by unity in coordinate 
space, but not in integrals, where it is retained in full form. The triplet density in the 
uniform limit approximation then reads 

P3(o, 1, 2, n) - i^p2(o, 1, n)p2{o, 2, n) ^ p'^pp {g{i, 2) - 1) , (2.21) 
p^ 

where g is the pair distribution of two ^He atoms, regardless of the position tq of the molecule. 

We can furthermore eliminate 5p2(0, 1, Q) by (5pi(0, Q) and 6u2, using eq. ()2.18p . and then 
make use of the 1-body equation ()2.16|1 . in order to arrive at 

= B{de^-P^^){deP2{0A,Q)) + B-y^{d^^ (2.22) 
r ^ 

+ p^p / d2 S(l,2) BL'^6u2(0,2,n) + p^ p — V?(5m2(0, 1, fi) 
J 2m 

+ |^(Vo^^^)(VoP2(0,l,fi))-p^p|^ / d2 S{l,2)VlSu2{0,2,n) 
+ ihp^p J d2 5(l,2)5u2(0,2,fi) . 

where the terms involving Vo stem from 6Cq^^^ /6u2{0, 1, Q). Note that the explicit reference 
to the external field V^'^^^\0,^l) has now been eliminated. This means that the 2-body 
correlation fluctuations are only driven by the 1-body correlation fluctuations, which in turn 
are the response to V^'-^^*''(0, ^l). In the above equation, S is the static structure function of 
^He in coordinate space, 

5(|ri - r2|) = 5(ri - r2) + p {g{\r, - r^l) - 1) , (2.23) 
S{k) = f d^re'^^Sir) . (2.24) 



The equations ()2.15|) and ()2.22j) can now be solved by expansion in plane waves and 
spherical harmonics. We define 

5pi(ro,fi) = E / 7^ e^'^'^YjMin) 5pjM{q) (2.25) 

V^(^^*)(ro,fi) = Y.j-^, e^'^'YjMin) vjlf (g) (2.26) 



J,M' 



(2vr) 
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6u2{ro,r„n) = ( 7^7^ e^'^'-Oe^P(r°-^)r,„(fi) a,„(k,p) (2.27) 

(27r)3f27r)3 



(27r)3 (27 

-^^],{kr)g,{k) ,(2.28) 



where r = ro — ri and cos a = v ■ Vt. 

If we restrict ourselves to an external perturbation potential that couples only to the 
rotational degree of freedom, i.e. V^^^\q = 0) = V^^^\ then translational motion is not 
directly excited. Since in CBF theory the molecule+helium system is regarded as being in 
its ground state before excitation, we are therefore calculating only the purely rotationally 
excited states, i.e. Spe^mil = 0) = 5p^,m. 

With the above transformations and after transforming from time to frequency, the 1- 
body and 2-body response equations ()2.15j) and ()2.22|1 become coupled algebraic equations. 
Eq. ()2.15|1 becomes 

huj5pjM{uj) + 2ViP{u;) = BJ{J+l)5pjM+ (2.29) 

+Anp^pB j:{-^Y'J Y;^,{n.,) ,Hp)«.™(o,p;-)(-i)^^(_^^' 



where we note that the translational current fluctuation in eq. ()2.15|) vanishes. Eq. ()2.22|) 
becomes 

ft^ 2 ^2 2 

^* + 2^d{i) + W " "'-(O-P!'^) = (2-30) 



X,jj. 



pp'-m ^'^^ p^ S{p) 



Here the static structure factor S{p) is defined in eq. fl2.24|) . The energy expression on 
the left hand side of ea. (I2.30p contains the free linear rotor spectrum B£{i + 1) and the 



Bijl-Feynman spectrum 4J] of bulk ^He, 



e{p) 



2mS{p) 

The symbols in the angular brackets follow from angular integration of spherical harmonics 



nil 1712 ni'^ \mi 1712 1713 



L{ii,e2,h) = h{ii,i2,i3) (4(4 + 1) + 4(4 + 1) -4(4 + 1))' (2.31) 

(24 + l)(24 + l)(24 + i) /4 4 4\' 



i:(4,4,4) 



47r \0 
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The expressions in round brackets are Wigner's 3-j symbols |45[. After eliminating agmiO, p) 
in eq. ()2.29p by using eq. ()2.30p . we use the summation rules for the Wigner 3-j symbols 
to further simplify eq. (|2.29p . It turns out that most of the angular quantum number 
summations are trivial and that the SpjM coefficients do not mix. This leads to the simple 
formula 

(fiJ(J+l) + SjH-M'^PjmM = 2V}p{uj). (2.32) 

Hence we have found the linear response of the density, Spj^j, to a weak perturbation of 
symmetry (J, M), Vja?"*- The excitation energies of the system are obtained by setting the 
perturbation potential to zero and solving eq. (j2.32p . 

htu = BJ{J+l) + (2.33) 

This has to be solved self-consistently in order to obtain the excitation energy ujj for given 
J (or energies, if there exist more than one solution for given J). These solutions correspond 
to the energies Ej = huj of the coupled molecule-hehum system, in the usual spectroscopic 
notation with J the total angular momentum associated with the solvated molecule. Sj(co') 
is the self energy 

' 2J + iyJ (27r)3 S{p) Bi{i+l) + e{p) + h^py2M -hu' ^ ' ' 

We note that T,j{uj) does not depend on the quantum number M anymore, and that therefore 
we will not observe any M-splitting of the R(0) line of HCN in bulk ^He. This splitting 
has been found in spectra of HCN in ^He dro plet s as reported in Ref. 0, where it has 
been attributed to the finite size of the droplets |l3l|. Our result for the excitation energy, 
eq. (|2.33|) , indicates that all M- levels are degenerate when HCN is embedded in uniform bulk 
"^He. Note that infrared spectra of HCN solvated in ^He droplets in a strong electric field show 
clear evidence of M-splitting Q]. The lack of M-splitting in our calculation of HCN in the 
homogeneous environment of bulk "^He thus indicates that the inhomogeneous environment 
of finite ^He droplets may indeed be responsible for the observed spectral splitting. 

In general, the self energy will be complex. Strictly speaking, eq. (|2.33p cannot be solved 
self-consistently in that case, and we can only speak in terms of the response of the system to 
the perturbation [i.e. laser field). The imaginary part of Sj(ct;), which is the homogeneous 
line-width, i.e. the inverse life-time of the state J, results from the contour integration 
/ dp when the energy denominator has a zero. At such quantum numbers p and £, energy 
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conservation allows for a decay of state J having energy huj into a lower rotational excitation 
with energy Bl{[ + 1), while exciting a phonon of energy e(p) and translational motion of 
the molecule of energy h^p^ /2M (momentum conservation) such that hujj = Bi{i + 1) + 
e{p) + h^p^ /2M. For all other combinations of p and a decay would not conserve energy. 
These decay channels are closed. 

The quantity L ()2.3ip contains a Wigner 3-j symbol as well as a rotational kinetic energy 
factor (the expression in the round brackets on the Ihs of eg. 12.311) . Combined, they obviously 
lead to the selection rules 

L{^l, 4) = , if ii + e2 + h uneven 

I/(£i,£2,^3) = , if £1,^2,^3 cannot form a triangle (2.35) 
^(4,4,4) = 0, if 4 = or £2 = 0. 

It follows that Sj(co') = for J = 0, i.e. the self energy does not renormalize the ground 
state energy. The coupling of the spectra in the energy denominator is mediated by the 
anisotropic pair probability distribution g{r, cos a), decomposed into its Legendre expansion 
coefficients. Furthermore, the spherical expansion coefficient gi'=Q does not contribute to 
Sj(cu). 

We note that in the self energy part of the spectrum ()2.33|) . hujj couples to free rotor 
states, to the free translational states, and to the Bijl-Feynman spectrum of helium. Al- 
though we should not over-interpret the meaning of the individiual terms in Sj(u;), in an 
exact expression for the correction to the rotational energy in the helium environment we 
would expect to find a coupling to renormalized molecular rotations and translations. We 
would also expect to see coupling to the exact energy spectrum of ^He, instead of to the 
Bijl-Feynman spectrum. Going beyond the uniform limit approximation for the probability 
densities might improve the molecular rotation spectrum in these respects, as has been found 



for other excitations in helium. These include calculations of the bulk helium spectrum 24 1 

. Alternatively (and much easier), we 
can choose a phenomenological approach and try to improve the self energy of eq. ()2.33p by 
using any one or a combination of the following replacements in the energy denominator of 
the self energy: 

1. In the following, we will always use the experimentally measured excitation spectrum 
instead of the Bijl-Feynman spectrum. 
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2. We can use the dispersion of translational motion of HCN and DCN in bulk ^He, 
fi^p^ /2Mqs instead of the free dispersion. Since we don't know of any experimental 
value for M^s (which would be a tensor quantity), we use the bare mass M. 

3. We can use hue self-consistently instead of Bi{£ +1). In this case, we solve hu = 
B£{i + 1) + S£(co') for hu angular quantum number i, and the solution hui replaces 
Bi{i + 1) in the energy denominator of Sj(co') for the next iteration; the procedure is 
iterated until convergence is reached for all huj. In the case of the calculation of the 
effective mass of impurities in ^He, this phenomenological approach was shown to to 
improve agreement with experimental results js^. However, we will see that in case 
for molecule rotations in ^He, for given J we can have several solutions uj of eq. ()2.33p . 
see appendix We minimize the ambivalence associated with this procedure and will 
not use this phenomenological improvement of the self energy. 

We discuss the dependence of the results on these phenomenological "improvements" in 
appendix O where B^q/Bq is calculated for various combinations of replacements 1, 2, and 
3. For the rest of the paper, we apply only replacement 1. 

A related concern is the missing of decay channels where localized "^He excitations are 
generated instead of a bulk helium excitation e{p). Localized layer phonons and rotons 



have been calculated |47| and observed j48| for helium adsorbed to graphite sheets, and 
localized vibrations calculated for helium adsorption on aromatic molecules j4^. Since 
the rotation of a molecule in ^He involves a correlated motion of the molecule and the 
surrounding '^He atoms, it can be regarded as involving a localized "layer" excitation of the 
^He when observed from the molecule frame. One significant difference from layer excitations 
of helium adsorbed to an extended substrate is that here the molecule "substrate" is so light 
that its motion must be taken into account (the rotational motion has been seen to have 
an influence on the vibrational energies for ^He adsorbed on the benzene molecule 
However, decay into channels other than bulk helium excitations, such as the localized 
molecule-helium excitations themselves, is beyond the ansatz of eq. (j2.7p and the uniform 
limit approximation, as we have pointed out above. Deriving and solving the CBF equations 
in the frame of the molecule would allow coupling to localized excitations, although it would 
considerably complicate the CBF equations. An extension in this direction might allow 
analysis of the rotational dynamics of heavier molecules such as OCS and extraction of the 
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moment of inertia renormalization deriving from coupling to some adiabatically following 
helium . However for the light HCN molecule, our results below show good agreement with 
the experimental rotational constant in helium, indicating that such coupling to localized 
excitations is not important in this case. 



B. Linear Response and absorption spectrum 

Once we have derived the excitation energy of molecular rotations in ^He from a linear 
response approach, we can also obtain the dynamic response function xi^)^ which is the 
linear operator relating a weak perturbation V"*^*^^*-* of frequency to and the response of the 
probability density 6pi: 

x(u;) = x'{^) +'ix"{^) consists of a real part yjui) describing dispersion, and an imaginary 



part x"{^) = 'S'(co') describing absorption |42L l5nl| (note that our definition of S^u) differs 



by a factor of tt from the definition of Ref. l42l ). However, we cannot simply identify the 
inverse of the expression in the bracket in eq. ()2.32j) with because Sp{uj) is the Fourier 

transform of the complex density fluctuation 6pi{0,Q). The physical density response in 
linear order is given by the real part of 5pi(0, Q): 

6pi{0,n;t) = ($(t)|pi(0,fi)|$(t)) = ^e6pi{0,n;t). 

Here pi(ro, fi) is the molecule density operator, which is given in coordinate space by (5(ro — 
rQ)5(fi — Q'). The expectation value of pi(0,f2) is the probability to find a molecule at 
position fq and orientation Q. 

To obtain x{^) from eq. ()2.32|) . we first note that 

SpjMiu;) = ^[6pjMiuj) + i-iy'5p*j^_^{-iu)] 

where we used the fact that V^'^'^'^\il;t) is real. With relation ()2.32j) we find 

Spjuico) = [Gj{co) + G}{-u)]v}^'\co) , 
where Gj{uj) is the resolvent 

Gj{uj) = [BJ{J +1) + Ejiuj) - nuj]-\ (2.36) 
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Since Gj{uj) is real for u < 0, we obtain for the dynamic response function 

Xj{oo) = Gj{uj) + Gj{-uj) . (2.37) 

From this the absorption spectrum of a rigid hnear rotor exposed to dipole (J = 1), 
quadrupole ( J = 2) etc. radiation of frequency u can be obtained as 

Sj{uj) = '^mxjioo) = $jmGj(cu) . (2.38) 
C. Marriage of DMC and CBF 

Formulations of the ground state Euler-Lagrange equations ()2.5|) which are suitable for 
numerical solution have to take advantage of the symmetries of the system under considera- 
tion. In our case this means translational symmetry and rotational symmetry around the axis 
of the linear molecule. Unlike the corresponding CBF calculation of excitation ()2.14|) . the 
ground state equations ()2.5p cannot be linearized, due to the strongly repulsive interactions. 
Consequently both their formulation for a specific symmetry and their numerical solution, 
are more demanding than the calculation of excitations. Nevertheless, the calculation of the 
self energy T,j{uj) (j2.34|) does require knowledge of some ground state quantities, in particu- 
lar of the ^He-'^He and the "^He-molecule pair distribution functions g{l, 2) (eq. (j2.23p ) and 
of (7(0,1, fi) (eq.EM)- 

The ^He-^He pair distribution function g{l,2) is the Fourier transform of the static 
structure factor S{k). For bulk ^He this has been obtained with great accuracy from neutron 



scattering experiments 
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52l |. S(k)h.a.s also been calculated using hypernetted chain / 



Euler-Lagrange theory (HCN/EL) |32| and DMC |53|. We have used the S{k) at T = OK 



from Ref. as well as the experimentally determined S{k) at T = IK from Ref. 0. These 
give essentially identical results for the rotational excitation energies, i.e. the results are 
independent of the finer details of S{k). We note that doping ^He atoms with a single 
molecule will cause only a change of S{k) on the order of 0(1/A^). Therefore we can safely 
use the S{k) of pure "^He in the expression for the self energy (j2.34|) . 

We additionally need to calculate the ^He-molecule pair distribution function g{0, 
Here Diffusion Monte Carlo (DMC) is of use: DMC is easy to implement for the calcu- 
lation of ground state properties, and since it does not require prior specification of sym- 
metries, one DMC implementation can be applied to any molecule-^He system with only 
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little modification. Hence we shall employ DMC for calculation of the ground states instead 
of solving equations ()2.5|) . This effectively avoids the difficulties of solving the non-linear 
Euler- Lagrange equations in a ground state calculation. We therefore use CBF theory only 
for excited states. The combined procedure can be summarized as follows: 

Step 1: DMC for calculation of the "^He-molecule pair distribution g{0, 1,^1); 

Step 2: CBF for calculation of rotational excitations huj and the corresponding absorption 
spectrum, using as input the ^He-molecule pair distribution g{0, 1, Q) obtained in 
step 1 and the ^He-^He pair distribution taken from experimental neutron scattering 
data 

The energies huj reported for HCN and DCN in this work (section llV Ajl are obtained using 
these two steps. 

In addition to the approximate calculation of excitation energies and life-times, CBF 
provides us with calculation of an excitation operator 6U for which 5f/|$o) is a good ap- 
proximation of the excited state wave function. This raises a potentially useful option for 
further improvement of energy calculations in these systems by direct means. We note that 

the representation of an excited state in terms of an excitation operator that is made in 

fl 

CBF is conceptually similar to the representation made in the POITSE approach |12|. In 
POITSE, the excitation operator provide input for a zero temperature imaginary time cor- 
relation function calculation from which the corresponding excitation energy is obtained by 
inverse Laplace transformation. In CBF, the excitation operator is one of the outputs of 
the calculation, and it is normally discarded. Finding the appropriate excitation operator 
for a POITSE calculation can be a hard problem in some systems. Therefore, knowledge of 
a good excitation operator deriving from a high quality CBF calculation may help consid- 
erably in reducing the computational expense as well as in simplifying the inverse Laplace 
transformation of a POITSE calculation. By using equations ()2.13p . ()2.27|) . and ()2.30p one 
can show that within CBF the 1-body term of 6U in Eq. ()2.7|1 is proportional to 

6ui{n)r^YjM{n) (2.39) 

i.e. the free rotor wave function, corresponding to an excitation energy BJ{J+1). Thus, it is 
the 2-body terms of 6U which are responsible for the reduction in value of effective rotational 
constant i?efr below the free rotor value B. To date, POITSE and related calculations for 
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rotational excitations of molecules in helium clusters III 
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have used only 1-body 



excitation operators of the above form. We therefore propose that in future implementations 
of spectral evolution methods such as POITSE, one employ the CBF excitation operator 5U 
of Eq. (j2.7j) . In this situation, the output of CBF, 5U , may then be used as the input to a 
third calculation step, namely 

Step 3: CBF provides the excitation operator for a POITSE calculation of the exact exci- 
tation energies huoj. 

We expect that because of the incorporation of molecule-helium correlations into the excita- 
tion operator within an exact calculation methodology, this should provide an improvement 
over the present calculations that terminate after step 2. 



III. RESULTS: GROUND STATE 



The implementation of DMC for a single linear molecule surrounded by He follows 
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57l |. treating the molecule as a rigid body with both, rotational and translational degrees of 
freedom. The difference is that here the system is confined to a simulation box of appropriate 
size and periodic boundary conditions are applied. The simulation box moves with the 
molecule such that the latter is kept in the center of the box (but the box does not rotate 
with the molecule). The size s of the simulation box can be either adjusted such that («) the 
system consisting of 256 ^He atoms and a single HCN or DCN molecule are in equilibrium, 
i.e. that the ground state energy is minimized with respect to variations of the box size; 
or such that {it) the ^He density reaches the asymptotic equilibrium value p = 0.022A^^ 
furthest away from the molecule (the edge of the simulation box). In the first method, the 
calculated quantity (the total energy) changes quadratically with the change of s, and in 
the second method the calculated quantity (the asymptotic density) changes linearly. The 
first method is thus more susceptible to errors by construction. The uncertainty in the total 
energy is largely due to the cut-off of the "^He-^He interaction potential at large distances 
(see below). For this reason, we chose to adjust s by the second method. In order to avoid 
excessive amount of calculations to find the equilibrium density (zero pressure), we choose 
only 3 box sizes, s =22. 5A, 23. OA, and 23. 5A. We found that the s =23. OA simulation yields 
edge densities closest to the equilibrium bulk value. Using one of the other box sizes did 
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not change our results for the rotational excitation energies within the statistical error. We 
have used a time step of dt = 0.15 niK for the imaginary-time evolution to the ground state. 
We have doubled dt and again have obtained the same result for the rotational excitation 
energies, thereby verifying that the DMC energies are free of finite time step bias. 

Ground state expectation values {i.e. g{r, cos a)) have been calculated with pure estima- 
tors using descendant weighting of importance sampled DMC according to the approach of 



"^T = exp - 



N 

^n[^^(|ri - ro|,cosai) + ^u'^\\ri - rj\) 

i=l i<j 



Ref. |58|. The trial wave function used here for the importance sampled DMC has the form 

(3.1) 

with the molecule-^He correlation m^^'' and the ^He-^He correlation ^2^'' given by 

u^"^^ (r, cos a) = — (c/r)^ (3.2) 
4^)(r) = -{b/rf (3.3) 

with c = 7.392 A and b = 2.670 A. The precise from of the trial function is not important 
because we use descendant weighting for obtaining unbiased values for g{r, cos a). Such 
an isotropic trial function was found to be adequate for previous important sampled DMC 
calculations for HCN in small clusters [l^. This is expected from the weak anisotropy of 
the HCN-He interaction. Fig. ^ shows contours of the molecule- helium interaction potential 
Vx for HCN-^He. For computational efficiency we introduce a cut-off for both the "^He- 
^He interaction and its correlation u'^^ at a radius = SA, and replace u'^^ by a smooth 
function 



U2^\r) = uf\r) - uf\rc) - (r - - — - 



du>F\r 



(3.4) 



dr 

For completeness, we report the ground state energetics of HCN in bulk ^He obtained 
within these calculations. In order to correct the total potential energy for the error in- 
troduced by the cut-off, we assume a homogeneous "^He equilibrium density p at zero 
pressure, and approximate the missing contribution to the total ground state energy, 
EcoTT = (p/2) d^rV}ie{r). This is not a highly accurate correction, but the exact value 
of the ground state energy is immaterial for our calculation of g{r,cosa). For the equilib- 
rium density p = 0.22A ^, the DMC sampling yields an uncorrected ground state energy 
of E'/N = -7.35 ± 0.006K. The correction is E^^rr/N = -0.95K per ^He atom. Thus we 
find a total energy of approximately E/N = — 8.3K for both HCN and DCN. The chemical 
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potential of the molecule /x is the difference between the energy E of molecule and helium 
and the energy of pure helium, Eq = N x 7.2K at equilibrium. Hence we find ~ — 282K 
for HCN and DCN in bulk ^He. 

In Fig. 121 the pair distribution g{r, cos a) (eq. (j2.28|) ) is shown for HCN in bulk ^He, 
simulated by 256 ^He atoms in a box of 23. OA length on each side with periodic boundary 
conditions applied. For DCN, we used the same box size. The coordinates r and a are the 
radial and polar spherical coordinates in the center of mass frame of the HCN molecule, 
with the molecular symmetry the z axis. 

Due to the small anisotropy of the ^He-HCN and ^He-DCN potential, and the large 
zero-point rotational motion of the molecule, the pair distribution g{r, cos a) is only slightly 
anisotropic. In Fig.|21 we show the Legendre expansion coefficients ge{r) of g{r, cosa), whose 
Bessel transform is the quantity entering the calculation of the self energy ()2.34|) . In the limit 
of S ^ oo, the zero-point motion would completely delocalize the molecule orientation with 
respect to the surrounding ^He. In this situation, g{r, cos a) would be isotropic, g£>o{r) = 0, 
and therefore the self energy correction to i?efr would vanish, B^s = B. With the large but 
finite B value of HCN, the Legendre expansion coefficients gi>o{r) are not negligible. As 
can be seen from Fig. 01 the quadrupole coefficient g2{'r) is the main contribution to the 
anisotropy of (^(r, cosa) for HCN. 

We recall that for the derivation of the rotational self energy expression ()2.34|) . the uniform 
limit approximation was applied (see discussion in section lTl A|l . For the Legendre expansion, 
this translates into the coordinate space approximations gi=Q{r) ^ 1 and (7^>o ^ 1. While 
all higher expansion coefficients ge>o{r) never exceed values of 0.2 in absolute value, go{r) 
deviates from unity considerably, varying between and values of almost 2. However, 
since, due to the selection rules, is independent of ge=o, we see that the extent of 

angular modulations in the helium solvation density are consistent with the uniform limit and 
that this is therefore a good approximation for the purpose of calculating purely rotational 
excitations of a light rotor like HCN and DCN. We note that for the heavier linear rotor 
OCS. WMC a ..onge. .o.e .n.o.op,c .nte.ac^on .eUu. Q. .he ^.u,a. 
modulation in the first layer of helium around OCS is considerably larger than ^ for HCN 
(Fig. 121). Hence the expansion coefficients ge{r) will be considerably larger and use of the 
uniform limit approximation would be more questionable for OCS. 



22 



IV. RESULTS: EXCITED STATES 



A. Rotational Energies of HCN and DCN in ^He 

The excitation energies are obtained as the solutions ujj of eq. (j2.33|) . Unhke for a hnear 
molecule in the gas phase (where Sj(co') = 0), it is possible that more than one solution 
exists for a given J. In the next section, we will show that this this is actually the case 
for J = 2 and J = 3 (and presumably for higher J's). The existence of several solutions is 
not surprising considering that uoj are the approximate excitation energies of a mani/-body 
system. 

Table El lists the energies of the primary rotational excitation for J = 1,2,3. By "pri- 
mary" we refer to the excitation of lowest energy, when we find more than one solutions of 
eq. ()2.33|) . The occurrence of several lines for a given J is discussed in the next section and 
in appendix Also shown in table H] are the respective experimental excitation energies for 
HCN and DCN obtained by microwave spectroscopy 1]. Only the energy for J = 1 could 
be measured experimentally, because the helium cluster temperature of T =0.38 K is too 
low to allow appreciable population of rotationally excited states for this system. 

B. Absorption Spectra of HCN in ^He 

As discussed in more detail in section QV E[ the self energy Sj(co') is complex, and the 
excitations obtained from eq. ()2.33p are therefore not true eigenstates but decay as a result 
of the coupling to ^He excitation modes. This effect is observed in the molecule absorption 
spectrum Sj{uj), eq. ()2.38|) . in the weak field V^^^""*) of frequency u. 

In a spectroscopic experiment, the frequency of a microwave laser field is scanned to 
obtain the rotational spectrum. Since the wavelength is much longer than the size of the 
molecule, only the dipole component of V^'^^^^{il), corresponding to the J = 1 component, 
is non-negligible. As a zero-temperature method, DMC/CBF only describes excitations 
from the ground state to an excited state. Hence, with the dipole field v}'^^'^ acting on 
the molecule, we obtain only the J = ^ 1 rotational excitation(s). This corresponds to 
the -R(O) spectral line. Neither the J = 1 ^ 2, 2 — 3, . . . excitations corresponding to 
i?(l), -R(2) and -R(3) spectral lines, nor the de-excitations J = 1 ^ 0, 2 — 1, . . . correspond- 
ing to the P(1),P(2) spectral lines are obtained directly. However, one can go from the 
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ground state to J = 2, 3, . . . simply by directly applying perturbations Vjj^j , J = 2, 3, . . ., 
i.e. via quadrupole-, octopole-, etc. transitions. Unlike in experiment, it is much easier 
in our CBF calculation to apply these multipole perturbations rather than formulate and 
solve the problem in a finite temperature theory. It has the added benefit that the zero- 
temperature absorption spectra Sj{uj) resulting from application of dipole, quadrupole, etc., 
perturbations are simpler to interpret than finite temperature spectra, while still containing 
all the information about the energetics of the system. In Fig. HI we plot the resulting 
absorption spectra Sj{uj), J = 1,2,3, for HCN, where we have set both the field strength 
K/M*'* the dipole moment of HCN/DCN to unity (these factors will only scale the inten- 
sities of the spectra). The DCN spectra looks very similar. As we have pointed out above 
(section III A|l , we correct the energy denominator of the self energy in ()2.37|1 by using the 
experimental phonon-roton spectrum instead of the Bijl-Feynman spectrum. However, we 
have not further replaced Bi{i+1) by huji = Bi{i+1) + Tii{uj£). A detailed discussion about 
the effect of these and other phenomenological corrections can be found in appendix O 

Without "^He surrounding the molecule, we have Sj(u;) = 0, i.e. the spectrum is a delta 
function centered at the free rotor energy 

Sf''''\u) = [hu - BJ{J + 1) + ie]'^ = n6 {hu - BJ{J + 1)) . (4.1) 

In Fig. 01 the free rotor lines are indicated by dashed vertical lines. 

In the CBF approximation, the self energy Tjj{uj) is finite with both a real and an imag- 
inary part. The associated absorption spectrum Sj{uj) shows two kinds of features - sharp 
peaks and broader bands. We first analyze the sharp peaks. The origin of sharp peaks in 
Sj{(jj) are discussed in detail in appendix El We show there how an imaginary part that is 
small relative to BJ{J + 1) + 3?eSj leads to a Lorentzian peaked at the energy huj, that is 
obtained as the solution (or one of the solutions) of hu — BJ{J + 1) — 9fteSj(ci;) = 0. The 
energy huj can be associated with a rotational excitation of finite life-time, which decays 
into a combination of a molecular i < J state and an excitation of the helium environment. 
The width of the Lorentzian is given by QmTijlujj) (see section HV E|l . 

In Fig. m the spectra Sj{uj) show sharp peaks of increasing width and decreasing height 
with increasing J. This indicates that the coupling of the HCN rotation to the phonon-roton 
spectrum of bulk ^He is stronger for higher energies. The lowest molecular mode J = 1 has 
the weakest coupling, evidenced by the fact that Si{uj) is very close to the spectrum of a free 
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rotor at T = 0, i.e. it has a single sharp hne. In the next section we will obtain the effective 
rotational constant i?efr from this line and compare with the corresponding experimental 
measurement. The exact width of the J = 1 spectral line is subject to computational 
uncertainties related to the DMC ground state calculation, as explained in section HV El and 
therefore cannot be directly compared with the experimental line width. In contrast to the 
single peak seen for J = 1, the spectra for J = 2 and J = 3 show several distinct sharp 
peaks. As explained in the appendix El calculation of the position of a peak can result in 
several solutions. In some cases [2^, the associated peaks have very small weight, but for 
J = 2 and J = 3, we find two clearly discernible peaks. Detailed analysis of the origin 
of these two peaks is also presented below in appendix ^ The analysis shows that this 
two-peak structure of Sj{uj) is a direct consequence of the divergent density of states of "^He 
at the roton minimum and the maxon maximum. Coupling to these divergences effectively 
splits the single free peak into two, and shifts one of the resulting peaks below the roton 
minimum and the other above the maxon maximum. Both peaks lie very close in energy to 
the divergent density of states of the phonon-roton dispersion. Therefore the motion of the 
molecule can couple to many excitations and the molecule rotates in a dense cloud of virtual 
roton and maxon excitations. Because of energy conservation, excitation of real rotons and 
maxons is not allowed at the energies of the two peaks. If it were allowed, it would lead to 
immediate damping of the rotation and we would not see well-defined peaks. 

We consider now the origin of the broader bands of Sj{uj) in Fig. 0] These broader 
bands are seen as additional features in the spectra for J = 2 and J = 3, between the two 
peaks. This is more clearly seen in Fig. El where the absorption spectra are now plotted 
all on the same scale and are shown together with the density of states for the bulk-|-freely 
translating particle excitation spectrum e{p) + h'^p'^/2M (bottom panel). In the energy 
range E = 12.0 — 15. IK, S2 has a high energy wing that is clearly aligned with the energy 
corresponding to maxon-roton excitations in ^He plus a recoiling HCN molecule. This wing 
structure is thus a signature of efficient coupling of the molecule to high-energy excitations 
in ^He that lie between the roton minimum and the maxon maximum. Excitations of the 
molecule in this wing feature are virtual, i.e. the molecule remains in its ground state i = 0. 
Another roton-maxon wing results from the coupling of the rotons and maxons with the 
i = 1 state of the molecule. Hence this wing is shifted by Bi{i + 1) = 2B and corresponds 
to the generation of a high-energy ^He excitation together with translational recoil of the 



25 



molecule, plus a molecular rotational excitation i = 1, i.e. the molecule is now not only 
translated but is also excited rotationally to the i = 1 state. shows qualitatively the same 
features. In contrast, in the roton-maxon wings are negligibly small. The primary peak 
has almost all the strength of the spectrum, because the dipole field directly couples to the 
J = 1 excitation energy of the molecule, the energy of which is much lower than the roton. 
Thus, for J = 1 alone the absorption spectrum looks like a gas phase spectrum and can thus 
be described purely in terms of an effective rotational constant -Befr which determines the 
location of the single peak. 

In principle there is an infinite series of roton-maxon wings for each i, shifted by Bi{£+1), 
with decreasing strength. However, with increasing energy hu, multi-phonon processes pre- 
sumably become important. For example, in pure helium, these processes become important 
for energies above approximately 25 K, above which the dynamic structure function S{k,uj) 
is dominated by multi-phonon excitations. |6l| In our implementation of CBF theory, only 
one-phonon processes are taken into account. 

Two technical details of our calculations are presented in the appendices. The first is the 
necessity to introduce a cut-off in the Legendre expansion of g{r, cos a) (appendix EI)- The 
second is a comparison of the effects of making the various phenomenological corrections to 
the self energy discussed in section ITTXl fappendix IH ]) . 



C. Effective rotational constant B^s of HCN and DCN in ^He 

From the position of the single peak in the absorption spectrum 5*1 we can obtain the 
rotational excitation energy of the J = 1 excitation, hui, from which we can obtain an 
effective rotational constant of -Befr assuming a free rotor spectrum: 

This is the direct analog of the procedure used to obtain an experimental measurement of 
Bcs in Refs. Table |n] compares the effective rotational constant of HCN and DCN, 

obtained from J = 1 only in this manner, with the corresponding measured values of Ref. ^ 
The statistical error of B^g shown in the table is propagated from the DMC ground state 
calculation of g{r, cos a). The values of -Befr are in overall good agreement with the exper- 
imental values, agreeing to within 5% for both molecules, although the error bars (~ 2%) 
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are unfortunately too large to confirm the experimental determination of a slightly smaller 
(~ 1.5%) ratio B^s/B for the lighter HCN than for DCN. For both HON and DCN, the 
calculated values of -Bgfr are slightly larger than the experimental values. Such behavior of 
CBF theory to produce somewhat higher excitation energies than the corresponding experi- 
mental (exact) values has been observed in other cases j32]. One remedy for this is to apply 
phenomenological correction to all terms in the energy denominator of Tjj{uj) as we have 
explained above. The values given in table |n] were obtained by making such a correction 
only to the Bijl-Feynman spectrum for bulk helium, i.e., replacing this by the experimental 
collective excitation spectrum, but not modifying B^{^ + 1) or the free particle dispersion 
h'^p'^ /2M. The additional effect of these further corrections is summarized in table UTTl where 
we see a slight improvement of B^^ in its agreement with the experimental values is obtained 
by making a self-consistent replacement of Bl{[ + 1) by hu^. Further improvement could 
presumably be achieved by replacing the bare molecular mass M by the effective mass Meg 
of HCN and DCN moving in ^He, if these quantities were known. However, we have checked 
that realistic changes in these quantities would not change the qualitative behavior of any of 
our results. These checks and relevant details for implementation of the phenomenological 
corrections are provided in appendix O 

Given that the value of B^s obtained here for HCN in bulk ^He is in good agreement 
with that measured in large droplets (A^ > 1000 helium atoms), it is interesting to compare 
a.o wi* ... co„.espondrns values calculated fo. ..all C^cQ (no experimental mea- 
surements have been made yet on small clusters). As noted in section HI calculations of the 
J = 1 excitation by the POITSE methodology show that the resulting fitted value i?efr does 
not converge to the large droplet value by = 25, in contrast to the behavior of the heavier 
molecules such as OCS and SFr. For these molecules, Bp^ converges to the corresponding 
depict value bcfo. the .olvation sKell is co.p.eteQ y H Q. The.e a.e 

several possible reasons for this difference. First, our analysis for HCN in bulk ^He shows 
that a light rotating molecule leads to generation of collective excitations that are extended 
in space (phonons and rotons) instead of to formation of a local non-superfluid density in 
the first solvation shell that can adiabatically follow the molecular rotation[4]. The cluster 
size dependence for these two different mechanisms might reasonably be expected to be very 
different, with the coupling to extended modes requiring more than a single solvation shell 
to approach its bulk character. A second possible explanation is that the projection operator 
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A used in Ref. Q accesses a higher J = 1 state than the state associated with the rotation of 
the molecule, thereby causing the targeted excitation to overshoot not only the large droplet 
value but also the gas phase rotational energy. This effect was already seen in the POITSE 
excitation spectra of the smallest clusters calculated in Ref. Q, where multiple peaks were 
found, one of which was consistently above the gas phase rotational energy. For = 1, 
comparison with the corresponding result obtained by the collocation method js^ confirmed 
that this excitation is indeed a higher lying J = 1 level. The POITSE method relies on 
having a good projection operator A'^q which has sufficient overlap with the desired excited 
state. In this case one seeks a rotation of the molecule, but one that nevertheless involves 
considerable correlation of the molecule with the helium, as is evident from the second term 
in the CBF excitation operator, eq. ()2.7j) . In contrast, the projector that has been used 



in both POITSE 



IJ, |l7[ and related methods |l8l| to date is a free molecular rotor function. 



which corresponds only to the first term in eq. ()2.7p . This suggests that it will be worthwhile 
to use the full stationary CBF excitation operator A = 6U in a POITSE calculation, as we 
have already proposed in section III CI above. 

The situation seems to be opposite for heavier rotors in ^He, namely here the direct 
approaches by evaluation and inversion of imaginary time correlation functions can provide 
a better description of the rotational dynamics. Thus, POITSE and related approaches are 
able to obtain accurate values for B^s for OCS but our current implementation of 

CBF is not expected to be reliable in this case, as we have noted in section Hi Al We expect 
that working in the frame of the molecule would improve the CBF description to account 
also for the adiabatic following of ^He around such heavy rotors. 



D. Effective distortion constant D^s of HCN and DCN in ^He 

We can use sharp spectral peaks for higher J values and fit to the spectroscopic energy 
levels for a non-rigid linear rotor, BJ{J + 1)-D{J{J+1))^ Q, 

where D is the centrifugal 

distortion constant. However, this fit should be used with considerable caution, for two 
reasons. First, as we have seen in section IIVBI the deviations of the J = —>■ 2 and 
— > 3 transitions from an effective linear rotor are very large and have nothing to do with 
a true centrifugal distortion that might arise from a coupling of the molecular rotation to 
either a molecular or a localized helium vibration. In particular, these higher transitions 



28 



are split into two peaks which cannot both be fit by a simple phenomenological centrifugal 
distortion term. Second, both the statistical errors from the DMC and the errors stemming 
from the approximations used in CBF (see section IH] above) grow with J, leading to larger 
overall errors associated with the peaks for higher J values. Nevertheless, by direct analogy 
again with the experimental procedure of fitting to spectral line positions, we can make a 
phenomenological fit to the J = 1 spectral line together with the lowest spectral peak for 
each of the J = 2 and J = 3 spectra. The result is that because of the downward shift 
of the J = 2 and J = 3 lines induced by coupling to the roton-maxon excitations, we find 
a very large value of the fitted effective centrifugal distortion constant Des- Thus, e.g., 
for HCN, while we obtain a value B^s = uJi/2 = 1.266 cm~^ (1.822 K) from the J = 1 
line only, fitting the J = 1 and first J = 2 peaks yields B^s = 1.346 cm~^ (1.937 K) and 
Des = 0.040 cm~^ (0.058 K), and further fitting the J = 1 and first J = 2 and J = 3 peaks all 
together yields B^s = 1-320 cm~^ (1.899 K) and D^s = 0.035 cm'^ (0.050 K). These values 
of Des are vastly enhanced over the gas phase value of the centrifugal distortion constant 
for HCN, D = 2.9 x 10^^ cm^^ [65|, showing an increase of four orders of magnitude. 
Similar enhancements of several orders of magnitude have been observed in experimental 
fitted values of D^s to multiple spectral lines for heavier molecules 13] and no theoretical 
explanation for these large enhancements has been given. 

The fitting constant D^s for HCN in helium and the gas phase centrifugal distortion 
constant D measure different physical effects. D is the usual measure of distortion of the 
linear rotor spectrum due to the non-rigid nature of HCN, which results in centrifugal forces 
acting on the component atoms as the molecule rotates, and hence in increased moments 
of inertia and lower rotational energy levels. D^s is a measure of the deviation from the 
effective linear rotor spectrum caused instead, in the case of the weakly anisotropic HCN 
molecule, by the "back-flow" of the surrounding helium. As our CBF results clearly show, 
the "back-flow" effect on the energy spectrum is much bigger than centrifugal distortion of 
the bare molecule. Thus the observed enhancement factor with respect to the gas phase 
value, Des/D, is not signiflcant, and D can be neglected in the discussion of rotational 
spectra of molecules in helium. 

In the flrst column of table HVl we compare the effective distortion constant D^s of HCN 
(from fltting to our CBF results) and OCS (fltted to experimental data) in helium. We also 
show values of Z^cfr obtained for HCN in a pseudo-hydrodynamic limit model, discussed in 



29 



section flV Fl In the second and third column, we show the respective ratios D^s/B and 
Dcs/Bes, i-e. we normahze Des such that all (free) linear rotor spectra would collapse on 
the same curve J{J + 1). 

We see from table HVI that regardless whether we use either -Defr? or D^s/B, or Dcs/B^s 
as a measure of distortion of the linear rotor spectrum, the HCN spectrum of excitation 
energies deviates considerably more from the linear rotor spectrum than does the spectrum 
of the heavier OCS molecule. We also see that the ratio Dcs/Bcs is similar for OCS and 
HCN in the pseudo-hydrodynamic model. A possible explanation for this last observation 
is given in section HV Fl below. 

The large value of Dgfr for HCN calculated here is a direct consequence of the high density 
of states in bulk helium near the roton-minimum and maxon-maximum, which are missing 
in the hydrodynamic limit model. This high density of states gives rise to a downward shift 
of the lower component of the split peaks for higher J states, as discussed in section IIVBI 
above and explained in detail in appendix El It also explains the greater distortion of the 
linear rotor spectrum compared to the distortion measured for OCS, since the considerably 
lower energy rotational excited states of OCS do not couple as effectively to the roton and 
maxon states as the rotational states of the lighter HCN molecule. 



E. Rotation life-times and homogeneous linewidth 

The self energy Sj has a small, but finite imaginary part, which leads to a finite life- 
time T = l/QmTij of the rotational excitation, i.e. to homogeneous line broadening of the 
rotational absorption spectra. 53mSj results from the principal value integration / dp in 
eq. ()2.34j) that is made when the energy denominator vanishes at some momentum p = po'. 

c^^y 2i?V pI Ee'L{J,i',i)gUPo) 

' 2J + 1 ^ Sipo) de{po)/dp + h'po/M' ^ ■ ' 

with Po defined by 

B£{e + 1) + e{po) = huj . (4.3) 

The life-time is obtained by summing the contribution from all poles po- For calculation 
of the J = 1 excitation, Eq. ()4.3|) only has a solution for i = and hence there is only one 
pole, because hur < Bi{i+ 1). From the selection rules ()2.36|) . we find that £' = 1 and hence 
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we obtain the estimate of line width 



3 ^S{po) de{po)/dp + h^po/M 



53ms,=i = ;b'p^ ... \; ,; , .2 ... • (4-4) 



Unfortunately the value of the momentum at the poles for HCN and DCN is very small: 
Po = 0.19A~^ for HCN and po = O.ISA"^ for DCN. These momentum values are too small 
for the corresponding Legendre component of the pair correlation function gi{po) to be a 
reliable estimate. This can be seen by considering for simplicity the Fourier transform of a 
periodic function with period s. The corresponding wave number p is discrete with smallest 
non-zero wave number equal to Pmin = 2tc/s. In our case, s = 23. oA is the simulation box 
length. This results in a minimum wave number pmin = 0.27A~^, which is larger than the 
desired pole momentum values Pq for HCN and DCN given above. Furthermore, we have 
the limiting value gi{p) = 0{p) for p ^ 0. Hence the Legendre component at the pole giipo) 
is small, resulting also in a small value of 53mSj=i. This explains the small width of the 
J = 1 line evident in Figs. |^ and El However, since the statistical error of gi{p) for small p 
is of the order of gi{p) itself, we are not able to extract a reliable quantitative estimate of 
the J = 1 life-time and the associated linewidth. 



F. Hydrodynamic Limit 

Hydrodynamical models have been used to describe rotations of heavy molecules with 



large moments of inertia solvated in He j4 



20|. As we have noted already in the intro- 



duction, these models fail for light rotors like HCN when based on assumptions of adiabatic 
following. In this section we show that independently of any assumption of adiabatic follow- 
ing, any analysis of light rotor rotation involving hydrodynamic coupling to long wavelength 
helium modes cannot provide an adequate description of the coupled molecule-helium exci- 
tations dynamics because of the absence of coupling of rotational levels to ^He excitations 
of higher energy than the long wavelength phonon modes, i.e. to rotons and maxons. 

We can simulate a hydrodynamic description of the ^He environment by replacing p- 
dependent ^He quantities by their low-p expansion: 

e{p) hep. (4.6) 
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For simplicity, we keep the molecule-^He pair distribution g{r, cos a) we have obtained from 
the quantum-mechanical DMC calculation. Therefore, our toy model is not a true hydro- 
dynamical model, which would require calculation of g{r, cos a) for HCN solvated in a hy- 
drodynamic environment. It should be noted that this "pseudo-hydrodynamic" model does 
not assume adiabatic following of the ^He ^,^|. In the present context, "hydrodynamic" 
refers simply to the coupling to bulk helium modes with long wavelength. 

In Fig. ini we show the absorption spectra Sj, J = 1,2,3 for HCN that are obtained 
with this pseudo-hydrodynamic model. These spectra show only sharp peaks and no broad 
bands, as expected from the discussion in section IIVBI that assigned the broad bands to 
coupling to collective excitations in the maxon-roton region. Also, only a single spectral 
line is found for all three J levels. According to the analysis in section IIV Bl and in the 
appendix El (see also Fig. |Hlin this appendix), this is also consistent with the lack of coupling 
to maxon-roton states. Fitting the positions of the three spectral peaks results in a good 
fit to a linear rotor spectrum, yielding effective spectroscopic constants B^s = 1.376cm^^ 
(1.980 K) and Des = 0.00568cm~^ (0.00817 K), respectively, and a corresponding ratio 
value Bcs/Bq = 0.931. The reduction in rotational constant is significantly less than the 
experimentally observed reduction of 0.815, amounting to only ~36% of the experimental 
reduction. This large discrepancy with the observed change of the rotational constant, in 
contrast to the good agreement achieved in section IIV CI from coupling to the true helium 
excitation spectrum further confirms that the maxon-roton regime of excitations has a strong 
influence on the rotation of light rotors in ^He. It is therefore essential for a reliable model 
to either directly (as in the present CBF approach) or indirectly (e.g., as in the POITSE 
approach) allow for coupling of the rotational levels with those excitations. 

The effective distortion constant Dcs is significantly reduced with respect to the full 
CBF result (see table HV]). while its ratio Dcs/Bcs is similar to the corresponding ratio 
for OCS. This suggests that the effective distortion constant of OCS may be rationalized 
as resulting from coupling to long-wavelength phonons. Within the present CBF analysis, 
this is consistent with the observation that for OCS with a gas phase rotation constant 
B = 0.2029cm^^ (0.2920 K), the roton excitations are too high in energy to couple effectively 
to the molecular rotation. However, for OCS the local coupling to helium^] needs also to 
be taken into account for a full analysis, as discussed in section Hi Al above. Consequently a 
consistent analysis of both B(,g and D^s for OCS will most likely benefit from reformulating 
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the CBF theory in the molecular frame, as suggested in section Hi Al 

Finally, it is interesting to note that this pseudo-hydrodynamic model severely underes- 
timates the change of the rotational constant for HCN, as opposed to the over-estimation 
for HCN that was obtained from the previous hydrodynamic model of Ref. Q that assumed 
complete adiabatic following of the molecular rotation by helium. 

V. CONCLUSIONS 

In this paper, we have derived the dynamic equations for molecular rotations in bulk ^He 
within the formalism of Correlated Basis Function (CBF) theory and applied them to HCN 
and DCN in superfiuid bulk ^He at T = 0. For that purpose we have combined Diffusion 
Monte Carlo (DMC) calculations for the required ground state properties with the CBF 
theory for excitations. Energy levels, absorption spectra, and spectroscopic constants for 
rotational excitations of the HCN and DCN molecules were calculated from this combined 
theoretical approach which allows for coupling to collective ^He excitations. Our results 
for the effective rotational constants of HCN and DCN are seen to be in good agreement 
with the corresponding experimentally determined values The CBF values are slightly 
higher than the experimental results (by ~ 5% of B, corresponding to ~ 25% of the re- 
duction B — -Befr); with about half of the difference being contained within the statistical 
error. We saw that the CBF values could be improved with systematic incorporation of 
more phenomenological input to the self energy. The statistical error derives from the DMC 
calculation of the molecule-helium correlation function and is hard to reduce further without 
imposing excessive computational requirements (the required sampling grows as the inverse 
square of the statistical error). For the present calculations with HCN and DCN, the sta- 
tistical error of DMC is unfortunately too large to determine whether the experimentally 
observed small isotope effect (~ 1.5%) is correctly predicted by CBF. 

An attractive feature of the CBF approach is the ability to calculate the full microwave 
absorption spectra at zero temperature. We calculated the spectra of the dipole, quadrupole, 
and octopole transitions of HCN, corresponding to J = ^ 1, ^ 2, and 0^3. The 
J = ^ 1 transition is found to be very sharp and the dipole spectrum seen to have 
almost no features apart from the single Lorentzian peak centered at the J = 1 excitation 
energy. In contrast, the J = — ^> 2 and J = ^ 3 transitions show weak phono n-maxon- 
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roton bands as well as secondary peaks. Both of these features are seen to be caused by 
the strong coupling of the molecule rotation to the roton and maxon excitations of "^He. 
This strong coupling was also seen to be responsible for the large values of the effective 
distortion constant D^s that result from fitting the primary peaks of the rotational excitation 
spectrum to the effective non-rigid linear rotor energy level expression BJ{J+ 1) — D{J{J + 
1))^. The importance of the phonon-maxon- roton spectrum was further highlighted by a 
comparative calculation where the rotational excitations are calculated with coupling to a 
phonon dispersion mode alone f section |IVF|) . In this pseudo-hydrodynamical model that 
lacks roton and maxon excitations, a much simpler absorption spectrum Sj was found 
that possesses only a single peak for all J values and no broad side-bands. The resulting 
empirically fit rotation constant is considerably higher than the experimental values, and 
smaller values of the effective distortion constant are seen. This shows that the coupling 
to the roton and maxon excitations of helium increases the deviation from the linear rotor 
spectrum. This coupling is strong for HCN and other light molecules, due to the vicinity of 
the J = 2 and J = 3 rotational levels to the roton energy. 

A key feature of our CBF results is their demonstration that the coupling to phonons and 
rotons of the bulk helium environment accounts quantitatively for the observed reduction of 
the effective rotational constant Bcs for HCN and DCN. For these light molecules, coupling 
to localized modes, e.g., as manifested by adiabatic following of some fraction of the first 
solvation shell helium density |^, can therefore be at most a very minor effect in the reduction 

of BeS- 

As far as methodology is concerned, the combination of DMC and CBF employed here fa- 
cilitates the calculation of excitation energies in the CBF approximation because of the ease 
of implementing DMC, which provides the ground state pair distribution 5f(r, cosa) needed 
for a CBF calculation. On the other hand, DMC is computationally much more expensive 
than the alternative of a full CBF calculation of both ground state and excitations. Even 
after extensive sampling, the statistical error of the ^He-HCN and ^He-DCN pair distribu- 
tion function g{r, cos a) was too large to be able to detect a statistically reliable difference 
between the rotational constants of the two isotopic species. Furthermore, with currently 
feasible simulation box sizes containing = 256 ^He atoms, we cannot reliably account for 
the long-range (small wave-length) correlations which are needed for the calculation of the 
homogeneous line width (section llV E|l . Ground state CBF, in contrast, is particularly reli- 
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able for long-range properties, while it yields only approximations to short-range properties, 
like the peak density of the first shell of ^He around the correlation hole of a ^He atom or a 
molecule. Thus, neither of these two approaches alone provides all the required ingredients 
to accurately obtain both the very small isotopic dependences of rotational constants and 
line widths for the HCN and DCN isotopomers. 

In this work, we have considered the simplest implementation of the CBF analysis for 
molecular excitations in ^He, by assuming an infinite bulk ^He matrix. However, matrix 
isolation spectroscopy experiments are performed in droplets consisting of a few thousands 
of "^He atoms. As noted above, the inhomogeneous environment has been shown to 
possibly cause inhomogeneous line broadening and may be responsible for the M-splitting of 
the observed R(0) line for HCN jj]. Another simplification made in the current CBF analysis 
was the disregard of coupling of rotation and translations of the molecule. This is justified 
at T = 0. At finite temperatures however, translational excitations will be populated, 
and these provide another source for inhomogeneous line broadening. The present CBF 
calculations can be generalized to molecules embedded in a finite quantum cluster which 
would allow quantification of the effect of a long range inhomogeneous helium environment 
(as distinct from the inhomogeneity in the local solvation shell around the molecule which 
is incorporated in this work) on the rotational dynamics of molecules. Finally we note that 
extension of the CBF approach presented here to heavy rotors like OCS and SFg may be 
feasible if the minimization of the action integral eq. ()2.8|1 is performed in a frame rotating 
with the molecule, thus allowing also for coupling to ^He excitations localized around the 
molecule and for adiabatic following of some local ^He density. 
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APPENDIX A: SELF-CONSISTENT SOLUTION OF EQ. (p35|> 



The absorption spectra Sj{lj) of eq. ()2.38p and ()2.36p are complete descriptions of the 
spectrum in that they contain all the information about excitation energies that CBF theory 
can provide. However, it is instructive to take a closer look at precisely how the sharp peaks 
in Sj{u!) arise from ea. (I2.38|) and ()2.36p . The present discussion follows closely the discussion 
of appendix A in Ref. l2a 

With the abbreviation 

7j(cj) = BJ{J + 1) + sReEj(cj) 

we can write the spectrum as 

c ( \ ^ '^m^jj^) 

Hence sharp peaks, i.e. long life-times of excitations 0| of energy Huq, occur when 77(0^0) = 
huo and QmT,j{u) is small. In this situation, Sj{u) is small everywhere except near u = ujq. 
In the vicinity of this region we can expand 



7j(a;) — Tiu} = ajh{uj — luq) with aj 
and obtain a Lorentzian centered at uq: 



SAcu) 



dhuj 



- 1 



UJ=UJQ 



h'aAoo-uJor + ey 

Here we have assumed that Qm'Ej{uj) varies very little in the region close to Uq and can con- 
sequently be replaced by ej = 53mSj (tuo). The weight of the peak is obtained by integration 
of the peak 

71 

dhu Sj{uj) = r- 



\aj\ 

Hence we find the position of a peak by solving the equation 

■jj{uj) — hu = 

for one or several roots cjj. We obtain the width of the peak from QmJ]j{uJi) and its weight 
from of \d'yj{uji)/duj — All this applies only when Qm'Ej{uJi) is small. 

In Fig.lSlwe show 'yj{uj) — hu as function of hu for J = 2. In the range shown in the plots, 
'yj{uj) — huj has two roots, cu^^^ and uj^'^\ indicated by black points. Since the imaginary part 
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of is small, 5'j(ci;) has a sharp peak at uj^^\ see Fig. In contrast, '^mHj^uj^'^^) is 

much larger, resulting in a broad peak at iu^'^\ Furthermore, both the real and imaginary 
part of T,j(uj) vary significantly near uj^'^\ so that the peak is no longer Lorentzian. A 
broad band can be seen between the two peaks in Fig. This broad band stems from the 
large values for '^m'Ej{uj) when hu lies in the roton-maxon band of the density of states, 
which is large between the two extremas of e{p) + h^p^/2M (see lowest panel of Fig. Ej). 
It is easy to show that the density of states as well as 53mSj(co') diverges as the inverse 
square root of the energy at the extremas of e{p) + h^p^ /2M (see lower panel of Fig. IH)). 
The real component '^eLj{uj) also diverges as the inverse square root, but it does so on the 
"outer" sides of roton-maxon band (see upper panel of Fig. |H)). These two divergences of 
3?eSj(u;) and hence of 7j(co') — huj at the roton-maxon band are clearly responsible for the 
occurrence of two roots uj^^^ and uj^'^\ In the pseudo-hydrodynamical model presented in 
section HV Fl we retain only the linear phonon dispersion and there is no roton-maxon band. 
Consequently, 9fteSj(co') does not diverge anywhere and we find only a single peak for each 
J in the pseudo-hydrodynamic calculations. 

APPENDIX B: £ CUT-OFF FOR c/^(p) 

Since g{r, cos a), and therefore gi{p), is affected by statistical noise, the self energy Sj(a;) 
is also affected by this. For large £, geXp) is small and the noise will exceed the true value of 
giip). But Sj(co') is a functional of gj{p), i.e. for large the summation over £ in eq. ()2.34|l 
adds only noise to Sj instead of converging. For that reason we introduce a cutoff £cut to 
geip), such that g^ip) = for £ > 4ut- In Fig. [3 we show the ratio Bcs/Bq for HCN as 
a function of the cutoff £cut, where the phenomenological self energy with the experimental 
rather than the Bijl-Feynman excitation spectrum was used in the denominator. Fig. [7| 
shows clearly that for HCN the largest contribution to Sj(co') comes from i = 2, i.e. the 
quadrupole deviation from a spherical distribution around the molecule. Beyond i > 4, 
geip) contributes very little to Sj, therefore we choose 4ut = 6. The uncertainty associated 
with £cut is much smaller than the statistical error of Bss/Bq that is propagated from the 
error of g{r, cos a). 
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APPENDIX C: CORRECTION TO T.j{io) 



The self-energy ()2.34|) was obtained by allowing for fluctuations of 2-body correla- 

tions and using the uniform limit approximation. As mentioned above, we can try to improve 
without changing its analytic form, but instead by introducing a phenomenological 



energy denominator, obtained by i) using the experimental excitation spectrum |l5l . |5J| in- 
stead of the Bijl-Feynman spectrum, ii) using the effective mass of HCN or DCN instead of 
the bare mass in or iii) using hui instead of Bi{£+ 1), as well as by combinations 

of these corrections. 

For the last replacement, we have to solve eq. (j2.33p self-consistently not only for J = 1, 
but simultaneously for all J, because of the occurrence of hui in Sj. Hence we solve the set 
of equations 

hiuj = BJ{J + 1) + 3ieSj(u;j), J = 1, . . . , J^^ (CI) 

with 

y( ) p2 (4^)'/" ^ f dp L{J, f , i)gUp) .^.^ 

■'^ 2J+l^i (27r)3 S{p) hue + e{p) + ny /2M - huj ^ ' 

We take the real part of the self energy, assuming that the imaginary part is small, since only 
then we have well-defined excitations, albeit decaying ones. We note that for J > 1, each 
one of eqns. ()(^H) has more than one solution, but we restrict ourselves to the solution which 
we believe corresponds to the effective rotational excitation of the molecule, i. e. to the main 
peak of Sj. This correspondence can only be established up to J = 3, hence we have to 
restrict ourselves to J^ax = 3. For J > 3 we use BJ{J + 1). In our view, solving eqns. ()C1|) 
for all other solutions as well and retaining the imaginary part of T,j{uj) would stretch the 
validity of a phenomenological correction of T,j{uj) and is therefore not warranted. 

The resulting four combinations of different corrections of Sj(c(j)'s for HCN in bulk ^He 
are compared in table ITTTl Clearly, the replacement of the Bijl-Feynman spectrum by the 
experimental excitation spectrum constitutes a significant correction of Sj(co') and manages 
to reproduce the experimental values of i?cfr- On the other hand, the self-consistent replace- 
ment of Bi{i + 1) by huJe leads only to a minor further reduction of B^s, almost within the 
statistical error of T,j{uj). Hence, we do not apply the latter correction in our calculations. 

It is instructive to consider the effect of the molecular mass more carefully. Unfortunately, 
the effective mass M^s of HCN and DCN in ^He is unknown. CBF permits the calculation of 
effective masses but this would be beyond the scope of the paper. Therefore we used the bare 
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mass M in the denominator of eq. ()C2|1 for all our calculations. In principle, we can turn the 
argument around and compare the solution ujj of eq. ()C1|) for J = 1 with the value for uji from 
the experiments of Ref. Q. However, on the level of CBF theory implemented in this paper, 
a precise prediction of the effective mass cannot be made, because uji depends only weakly 
on Mofj. In Fig. El we show the ratios of calculated versus experimental effective rotational 
constant, Bcs/B^q, for HCN and DCN as a function of the effective mass ratio Mes/M, 
where Bes = 2uji has been obtained from eqns. ()C1|1 and ()C2|) . using the experimental "^He 
spectrum for t{p). The error bars in Fig. El are estimated from the statistical error of 5eff 
obtained from eqns. ()2.33|1 and ()2.H4|1 . with M^s = M. The curves intersect Bes/Bl^ = 1 at 
different values of M^^/M, implying widely different effective masses for the two isotopes. 
However, the error bars are very large, and for most of the range in Fig. El the dependence 
of -Beflf on Meff is not statistically significant. Note also that we have neglected the effect of 
coupling of translation and rotation on M^.^. Coupling would of course introduce directional 
dependence of M^s in the molecular coordinate frame. 
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TABLE I: Energies of the primary rotational excitation of HCN and DCN. CBF denotes the present 
calculations employing CBF theory for excitations combined with exact ground state quantitites 
calculated by DMC, and exp. refers to the experimental values of Ref. Q. 





HCN(CBF) 


HCN (exp.) 


DCN(CBF) 


DCN(exp.) 


J = 


1 


2.53 cm-i 


2.407 cm-i 


2.08 cm-i 


1.998 cm-i 


J = 


2 


6.64 cm~^ 




5.76 cm~^ 




J = 


3 


10.8 cm-i 




9.77 cm-i 
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TABLE II: Comparison of the calculated ratio Bcs/Bq of HCN and DCN with the corresponding 
experimental values Q]. CBF refers to the current CBF theory for excitations combined with exact 
ground state quantities calculated by DMC. Within the statistical error, the calculated ratios 
B^s/Bq for HCN and DCN cannot be distinguished from each other. The last column lists the 
corresponding moment of inertia increase (in uA^) in CBF theory and in experiment, respectively. 





CBF 


experiment [Ij 


AI (CBF) AI (exp.) 


HCN 


0.857 ± 0.019 


0.814 


1.90±0.29 2.61 


DCN 


0.863 ± 0.016 


0.830 


2.22±0.31 2.87 
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TABLE III: Comparison of the calculated ratio Bcs/Bq for HCN obtained with and without 
the phenomenological corrections explained in the text. The four entries correspond to the four 
possible combinations of corrections in the energy denominator of the self energy (i) the 

Bijl-Feynman spectrum (left column) or the experimental spectrum (right column), and (ii) the gas 
phase rotational energies (top row) or self-consistent solution of the rotational energies in helium, 
eqns. (fnT|) and (fn2|) (bottom row). 





Bijl-Feynman 


exp. spectrum 


Biil + 1) 


0.913 


0.857 




0.910 


0.841 
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TABLE IV: The effective distortion constant Dcs and the ratios Dcs/B and Dcs/Bcs calculated 
for HCN in helium by the present combination of CBF and DMC, compared to the corresponding 
ratios derived from experimental measurements for OCS in helium droplets Gg]. We also show the 



corresponding values calculated for HCN within the pseudo-hydrodynamical model of section HV Fl 
The last two columns give the gas phase reference values of the spectroscopic constants B and D 
for the two molecules. 





Des (cm~i) 


Des/B 


Des/ BcS 


B (cm-i) 


D (cm-i) 


HCN 


0.035 


0.0237 


0.0265 


1.478 \65] 


2.9x10-6 [65j 


OCS(exp. [66]) 


0.0004 


0.00197 


0.00546 


0.0732 


0.438 xlO"'^ 


HCN (pseudo-hydro. ) 


0.00568 


0.00384 


0.00412 
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-6 

2 4 6 

R[A] 

FIG. 1: Contour plot of the HCN-^He potential surface Vx{R, z) from Ref. [sj. Contour levels are 
shown at energy spacings of 5 K, with the outermost contour at -5 K, the next one at -10 K, etc. 
The linear HCN molecule is oriented along the z axis such that the hydrogen points in the positive 
z direction. R is the cylindrical polar radius. At z = 4.25 A and i? = oA, the potential attains its 
minimum value of -42.4 K. 
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FIG. 2: Pair distribution g{r, cos a) (eq. XTM ). between HON and ^UeN, for = 256. r is the 
distance between HON and a '^He atom, and a is the angle between the directional vector from 
HON to ^He and the HCN axis. The HCN-^He interaction potential (see Fig. ^ is defined such 
that hydrogen is located on the positive side of molecule axis, i.e. at a = 0. 
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FIG. 4: The absorption spectra 5j(u;), J = 1, 2, 3, for HCN in bulk ^He, where this is represented 
by = 256 ^He in a box subject to periodic boundary conditions (see text). The dashed hues 
indicate the corresponding rotational excitation energies of HCN in the gas phasej^^. The spectra 
have been broadened by a Lorentzian, by adding a small constant imaginary part of lOmK to the 
self energy T,j{uj). 
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FIG. 5: The spectra Sj{ul>), J = 1,2, 3, for HON in bulk "^He with all Sj{lo) now shown on same 
scale. For details of the representation of bulk ^He see Fig.^and text. The bottom panel shows in 
addition the dispersion curve e{p)+h'^p^ /2M and its density of states ("dos") [de{p) / dp+h'^p/Ai]~'^ . 
The vertical dashed and dotted lines indicate the onset of the roton-maxon band coupling with the 
£ = (dashed lines) and the 1 = 1 (dotted lines) rotational states of the molecule, respectively (see 
text). 
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FIG. 6: The absorption spectra Sj{uj), J = 1,2,3, for HON obtained with the "pseudo- 
hydrodynamic model" in which the bulk "^He is replaced by a hydrodynamic model fluid having 
only long wavelength (phonon) modes, i.e., possessing linear dispersion. The simulation is made 
with N = 256 ^He in a box subject to periodic boundary condition and a length of 23. oA. The 
dashed lines indicate the corresponding spectral positions for HON in the gas phasej^^. 
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FIG. 7: The effective rotational constant ratio B^s/Bq as a function of the cutoff £cut for HON 
in bulk ^He, approximated by = 256 ^He in a box subject to periodic boundary condition. We 
have truncated the Legendre expansion of the pair distribution g{r,a) at £cut- The expansion has 
reached convergence at a value slightly higher than the experimentally measured ratioQ] which is 
indicated by the horizontal line. 
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FIG. 8: Illustration of procedure for obtaining excitation energies from eq. (|2.3.S|) . The upper panel 
shows 7j(u;) — hto for J = 2, where Jjiuj) = BJ{J+ 1) + ^eT,j{uj). The points emphasize the zeros 
of 'yj{uj) — huj, which are the excitation energies for J = 2. For comparison, the dashed line shows 
7j(w) — hu), J = 2, from our pseudo-hydrodynamic model, see section HVFl where 'yj{u)) = huj has 
only one solution, leading to the single peak shown in Fig. |21 The lower panel shows 9mSj(6<j) 
for J = 2. Again, the dashed line indicates the corresponding result of the pseudo-hydrodynamic 
model. See text for a full discussion. 
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FIG. 9: The error of calculated versus experimental effective rotational constant, Bcs/B^s, for 
HCN and DCN as a function of the effective mass ratio M^s/M. That the intersections of the 
curves with Beg/B^^ = 1 occurs at the same values of Meg/M, is clearly not statistically significant 
because of the large error bars. The value of B^^ is taken from Ref. Q. See appendix [0 for a 
discussion of the phenomenological corrections used in this figure. 
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